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Abstract. 

In this paper we show how the formalism of Wu & Goldreich developed for 
weakly nonlinear white dwarf pulsations can be extended to the fully nonlinear 
ca se. Due to the co mpu tational simplification this affords over the simulations 
of lBrir.khill ll1992aD and Using KoesteH lEKVty . we are able to obtain the first 
ever fits of observed light curves. These fits allow us to obtain not only mode 
identifications but also to constrain the average depth of the convection zone 
and its sensitivity to changes in T e ff . We are thus able to probe the fundamental 
physics of convection in these objects. 



1. Astrophysical Context 

In the past decade, the seismology of pulsating stars in general ( "asteroseis- 
mology"), and of white dwarfs in particular, has produced remarkable results 
concerning the parameters and internal structure of stars. Besides determining 
very precise luminosities an d therefore distances to many of these stars ( "astero - 
seismological parallax", e.g. lBradlev h WingetJll994l : IWinget et al.lll994l . Il99lh . 
we have placed constraints on many aspects of their i nternal structure, su ch as 
their helium and hydrogen surfac e layer masses (e.g. iMetcalfe etaTl feoOOt and 
their core carbon/oxygen profiles f Metcalfdl2003l ) : in addition, we have provided 
empirical evidence that crystallization in cooling white dwarfs does in fact occur 
( Metcalfe et 111120041 : iMontgomerv fc Winged 119991 : ISalneterlll96lT ). Finally, we 
have recently made fundamental progress in understandin g how the observed 
eigen fr eq uencies sample the radial structure of our models ((Montgomery et alJ 
2003), and how the resulting "core/envelope symmetry" means that our models 
are actually very sensitive to the structure of the deep core. 

However, all of this progress in the field of white dwarf asteroseismology 
has involved the observation and modelling of their linear, sinusoidal pulsations. 
While this program of investigation has been very fruitful, it is somewhat ironic 
that about 2/3 of all white dwarf pulsators are large- amplitude and nonlinear. 
This sets them apart from most other classes of non-radially pulsating stars, 
which are usually small amplitude (e.g. 5 Scuti, (3 Cephei, 7 Doradus, and solar- 
like pulsating stars), and places them with the large-amplitude radial pulsators, 
such as the RR Lyrae and Cepheid stars. 

These observed nonlinearities contain information about the physical pro- 
cesses which produced them. Along with previous researchers, we believe these 
nonlinearities to be the result of th e interaction of the pulsations w ith the con- 
vection zone (e.g. iBrickhillll 1 992al lbL IWull200ll : Using fc Koester1l200lh . Physically, 
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the convection zone varies in depth and therefore heat capacity during a given 
pulsation cycle, with the result that it modulates the (assumed) sinusoidal flux 
flowing into it from below, producing a nonlinear flux variation at the photo- 
sphere. The amount by which the depth of the convection zone varies is inti- 
mately tied to the physics of convection in these stars, and therefore is a probe 
of this physics. Given that convection is even today somewhat unconstrained in 
the modelling of these and other stars (the "free" parameter a) , any progress we 
can make in providing independent observational constraints is of definite value. 

In the sections which follow, we show how a particular physical approxima- 
tion can be used to greatly simplify the description of the pulsation/convection 
interaction in these stars. Using models based on this, we apply this method to 
the observed pulsations of two different stars. These fits to observed light curves, 
and the subsequent derivation of convective and pulsation mode parameters, 
represent the first such results for white dwarf p ulsators, and provide important 
tests an d points of contact for both older (e.g. lB5hm-Vitense1ll958t ) and more 
modern ([Montgomery fc Kunkal[2004l : iKupkal llQQ^ theories of convection. 

2. The Model 
2.1. Assumptions 

While we are simulating the same basic physical system described in iBrickhilll 
( 1992aT ) , the set of assumptions which we adopt allows us to model and therefore 
fit light curves much more easily (see IWull200lh . The fundamental assumptions 
of this model are that: 

1. The flux perturbations beneath the convection zone are completely sinu- 
soidal. 

2. The convection zone is so thin that we may locally ignore the angular 
variation of the nonradial pulsations, i.e., we treat the pulsations locally 
as if they were radial. 

3. The convective turnover timescale is so short compared to the pulsation pe- 
riods that the convection zone can be taken to respond "instantaneously" . 

4. Due to the extreme sensitivity of the convection zone, only flux and tem- 
perature variations need to be considered, i.e., the large-scale fluid motions 
associated with the pulsations can be ignored. 

The third assumption is given weight by appeal to standard convection models 
of these stars, which indicate that the time taken for a convective fluid element 
to circulate from top to bottom of the convection zone should be less than a 
second. Since the pulsation periods of these stars are much longer than this, i.e., 
hundreds of seconds, this appears to be a safe assumption. While plausibility 
arguments for the other assumptions may similarly be made, it is perhaps better 
to see how well they allow us to model observed light curves before deciding on 
their viability ("the proof is in the pudding"). 
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2.2. Energetics 

We assume energy balance in the sense that the energy radiated away at the pho- 
tosphere (the top of the convection zone) is equal to the energy absorbed at its 
base minus the energy absorbed by the convection zone itself. Mathematically, 
we write 

F - F - d ® 

-fphot — ^base ^ j 

energy energy 

output mput absorbed 



i.e., the flux emitted at the photosphere equals the flux incident at the base of 
the convection zone minus the rate at which flux (energy) is absorbed by the 
convection zone. Physically, this can be thought of as the convection zone having 
a specific heat. For instance, in order for its temperature to be raised by 200 K, 
say, it must absorb a given amount of energy. 

Due to the assumption that the convection zone is instantaneously in quasi- 
static equilibrium, we can compute the term dQ/dt purely in terms of static 
models. As it turns out, this energy absorption rate depends only on the pho- 
tospheric flux (i.e., T e ff), leading to the following equation: 

-Fphot = F hase + tc — Jjp^s (1) 

where the new timescale tq = Tc{F phot ) describes the changing heat capacity 
of the convection zone as a function of the local photospheric flux. Thus, we 
have reduced the problem to a first-order equation in time. However, we must 
remember that although i^ase is assumed to be sinusoidal, it still does have a 
spatial dependence (-Fbase °c Re[e lujt Yi m (0 , <p)]), so equationHmust be solved on a 
grid of points across the visible surface of the star and the resulting fluxes added 
together to calculate the observed light curve. Fortunately, as a consequence of 
the simplifications we have introduced, this problem is still tractable. 

In Figure 1, we show a theoretical calculation of how this timescale tq is 
expected to vary as a function of T e g. For a given value of a, we see that tq 
changes by approximately a factor of 1000 as the temperature goes from 12000 K 
(the observed onset of WD pulsations) to 11000 K (the observed disappearance 
of pulsations). Since a typical pulsator may have excursions in temperature of 
several hundreds of degrees, we expect tq to vary greatly during the pulsations, 
leading to nonlinear light curves. 



3. Light Curve Fitting 

Using Figure a guide, we see that although tq is a strong function of tem- 
perature, its logarithm varies approximately linearly with temperature. Since 
5F/F = 4 5T/T, log tq is also an approximately linear function of F, and since 
we wish to make as few assumptions about convection as possible, we choose 

where tq and a are unknown coefficients and Fq is the equilibrium flux. 



4 



Montgomery 




1.1 1.15 1.2 1.25 

T eff /( 10000 K) 

Figure 1. Theoretical calculations of the thermal response timescale of the 
convection zone, Tn. as a function of T e g. The local mixing-length theory of 
Boh m fc Cassineilj(|1971|) has been used, for the indicated values of a. 

In addition to the convective parameters tq and a, we fit the amplitude and 
phase of the sinusoidal signal at the base of the convection zone (the period is 
assumed already known from the observations) as well as the inclination angle 
between the pulsation axis and our line-of-sight. 

3.1. The DBV PG1351+049 

In 1995 data were acquired on the DB V star PG1351+049 as p art of a Whole 
Earth Telescope (WET) campaign (e.g. lWinget et al.lll99ll . ll994h . The resulting 
light curve, which has been folded at the period of the mode, 489.335 s, is shown 
in the lower panel of Figure [3 where the crosses represent the data and the solid 
curve is our fit. We note that the fit is very good and reproduces the features in 
the light curve quite well. In the upper panel we show the convective timescale 
which we derive from these data (solid curve) as well as one computed using 
mixing-length theory (dashed curve). The fact that the slopes are similar is 
encouraging, since the theoretical curve can be moved up or down simply by 
tuning the value of a, while its slope is less sensitive to the choice of a. 

This star was observed again in May 2004 with the 2.1m telescope at Mc- 
Donald Observatory. This allows us to make another fit, which provides an 
important check since we would not expect the parameters characterizing the 
star itself to have changed, i.e., the average depth of the convection zone, to, and 
the inclination angle, should be the same, while the amplitude and phase of 
the mode might very well be different. We summarize these two fits in Table 1: 
we see that the amplitude of the mode has changed somewhat, but 9\ and To 
remain the same, i.e., we appear to be looking at the same star. 



Tablel. Derived parameters for PG1351+049 



year to (sec) 


Oi (deg) 


Amp 


I m 


1995 99.95 


59.4 


0.325 


1 


2004 100.83 


59.4 


0.268 


1 
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Figure 2. A fit to the observed pulse shape of the DBV star PG1351+049. 
3.2. The DAV G29-38 

The observations of the DAV G29-38, taken by S. Kleinman in 1988, are dis- 
played as the crosses in the lower panel of Figure El (note: the errors on this 
extremely high signal-to- noise data are actually much smaller than the crosses). 
Our fit to this data is given by the solid curve, which is nearly indistinguishable 
from the data except at the maxima and minima, where slight differences can 
be seen. In the upper panel, we show the derived dependence of tq as a function 
of T c g (solid curve) as well as a theoretical curve computed using mixing-length 
theory (dashed curve). We see that the slopes of these curves are nearly iden- 
tical over this region, again implying that we really are measuring something 
fundamental about this star, and about the physics of convection. 

4. Conclusions 

First, the fits we have obtained for the stars PG1351+049 and G29-38 are im- 
pressively good, implying that the physical assumptions which have gone into 
our models are realistic. Second, our results seem to make sense in terms of sim- 
ple (and naive) theories of convection, and they provide a means of calibrating 
such models. Third, the fits provide an independent method for determining the 
£ and m values of pulsation modes, something which is necessary in order to use 
asteroseismology to learn about the deep stellar interior. 

The ultimate goal of this research is to provide a detailed empirical map of 
how convection works in these stars, which in terms of their pulsations means 
mapping how tq varies as a function of T e g and log g. Since each pulsator es- 
sentially gives us the value of tq and its slope at a given T e g, by applying this 
method to stars with temperatures spanning the instability strips of the DAVs 
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Figure 3. A fit to the observed pulse shape of the DAV star G29-38. 

and DBVs we can derive this dependence empirically. In the future, this de- 
termination will provide an important contact for more sophisticated and phys- 
ical theories of convection which go beyond mixing-length theory (those with- 
out " free parameters" , e.g. iMontgomerv Kupkall2004l : iKupka Montgomery! 
2002]) and whose applicability will be to stars throughout the HR diagram. 
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